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. . Contrary to recent theoretical prediction, we show that the magnetic 

flux of a vortex in SU{2) model of two-gap superconductor is quantized 
in units of 2Ti/g, not 'iii/g. For the U(l) version of this model, the flux 
is quantized in units of 2-Ka/g where < a < 1. The parameter a 
depends on the masses and concentrations of the Cooper pairs in the two 
condensates. 



Abstract 



PACS number(s): 74.20.De, 74.25.Qt, 74.25.Sv, 74.90.+n 



Ginzburg-Landau models with two flavors of Cooper pairs have attracted 
increased attention in recent years [1-3] . The discovery [4] of the two-band su- 
^t^ ' perconductor MgB2 with surprisingly high transition temperature has sparked 

\^ , renewed interest in these models. Starting from a Lagrangian for a charged 

■^ ■ doublet order parameter exhibiting SU{2) symmetry, Cho [2] studied the flux 

>^ I quantization for a magnetic vortex. In his approach , the magnetic flux is re- 

lated to the boundary values of the gauge field at r = and ?■ = oo where 
r is the distance from the vortex center. The r = value is determined by 
j^ ' examining the equation of motion for the magnitude of the order parameter. 

Smothness of the order parameter at r = requires that this boundary value 
is nonzero. Consequently, the flux is predicted in Ref. [2] to be quantized in 
^2 I the unit iir/g which is double of the flux carried by the well known Abrikosov 

vortex in one-gap superconductor [5]. 
fj , In the present paper, we consider the same SU(2) model of two-gap super- 

conductor and calculate the magnetic flux for the vortex in terms of the Berry 
connection that is due to the texture induced by the vortex in the charge space 
J^ ' [2, 3]. In this way, we find that the vortex has a total fiux quantized in the 

JH I same unit as the Abelian vortex in one-gap superconductor. Moreover, the 

present method allows us to consider the fiux quantization in a more realistic 
C/(l) model of the two-gap superconductor. In this model the flux appears to 
be quantized in the unit of 2'K/g times a parameter that can take any value in 
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the interval (0, 1) dependent on the masses and concentrations of the Cooper 
pairs. 

The Lagrangian density has the same form as that of Ref. [2] 

£ = _|2?,^|2 + ^2^t^ _ ^(^t^)2 _ 1 j,2^ (1) 

We consider a magnetic vortex in a (2 + l)-dimensional spacetime. For a 
time-independent configuration, </) is a function of two space variables (polar 
coordinates (r, </?)). The object D^4' is the covariant derivative 

D,4,^{d, + igA^)(t) (2) 

where A^ is the Abehan gauge field given by 

Tfi 

A, = -A(r)a,(^ (3) 

where ^ is the phase entering the vortex Ansatz (j){r, ip) parametrized by 
1 1 / cos ^^-^ exp(— i?7i(^) 



where /(r) decreases monotonically from /(O) = tt to /(oo) = 0. 

We note that expression (3) has a form dictated by requiring that the two 
terms in Eq. (2) cancel each other as r — ^ oo so that the first term in Eq. (1) 
yields a finite energy [6] . 

Varying £ with respect to A^^ yields the gauge invariant current 



Ji^ = *5 



{D,^)+^-^^iD,cj)) (5) 



-r m 

A„ = - — 



Using spinor (4), this equation takes the form 

>=5VM(^-+^.) (6) 

where A^ = —-^^d^^ is the Berry connection characterizing the bundle of 
eigenstates (4). With the spinor ^ defined in (4), we have 

\l + cosfir)\d,ip (7) 

Let us consider Eq. (6) for a large circle C of radius R in the xy-plane. 
For large r, the current density j^ decays exponentially owing to the Meissner 
effect. Also the density p{r) -^ po as r ^ oo. Thus, Eq. (6) implies that 
A,y{R) -^ —A^{R) as i? ^ oo. Integrating this relation over C, we obtain with 
the use of the Stokes theorem the total magnetic fiux passing through C 

^ — curlzA(fx = m A^dx^ ^ ^ f ^v^Xi, (8) 

J Jc Jc 

as R ^ oo. The right hand side (R. H. S.) of this equation can be evaluated 
using Eq. (7) and noting that /(oo) = 



~ TTi r 1 TTi 

MR)^-7r i + cos/(i?) d^if^ — d^ip 

as i? ^ oo. According to Eqs. (8) and (9), the flux is given by 

27rTO 



$ 



5 Jc 



d^(j)dxi, 



9 



(9) 



(10) 



This result contrasts with the flux ^2i™ obtained in Eq. (7) of Ref. [2]. Note 
that the gauge field A{r) at the origin r — Q does not appear in our derivation. 

To understand this disagreement, we examine the curl of the gauge field (3) 
paying attention to the singularity at r = 



curhA = e^"'^aA^ 



m 1 dA{r) m 
g r dr 



9 



A{r)e^'-d^d,^ 



(11) 



The second term on the R. H. S. of this equation is proportional to curlz — 
where 0q is the unit vector of the polar coordinates. This quantity vanishes 
except along the line r = 0. Applying the Stoke's theorem to a small circle 
around the origin, one obtains [7] 



e^^dfj^dv^p — curl 



Vo 



2TT5'^{r) 



(12) 



Using Eqs. (11) and (12), the magnetic fiux becomes 



$ = / curLAdx = 



2TTm 



27rm 
9 



dA(r 
dr 



-dr 



A{r)6^{r^d^x 



A{oo) - A(0) 



A(0) 



2'K'm 



A{oo) 



(13) 



We see that the A(0) term appearing in Eq. (7) of Ref. (2) is cancelled in 
Eq. (13) by the singular term stemming from Eq. (12). According to Eq. (5) 
of Ref. [2], we have A{oo) = 1. With this value, Eq. (13) predicts a fiux 2'iTm/g 
in agreement with Eq. (10). 

In the differential equation for p{r) as derived in Ref. [2], the quantity A{r) 
enters in the combination [A{r) — £2£l±i]2 _^ A^{{]) as r ^ 0. This is to be 
compared with the form [A{r) — 1]^ pertinent to the one-gap Ginzburg Landau 
superconductor. Consequently, the requirement of smoothness of p{r) at r = 
imposes the boundary condition Air) ^ — 1 as r ^ [2]. This contrast with 
the boundary condition A{{)) = for the one-gap case [6]. 

Ref. [2] proposes to remove this anomalous boundary condition by the gauge 
transformation 



A,. 



A' = A,. 



'si 



(14) 



Recalling Eq. (3), this implies A'{r) = [A{r) + 1] leading to new boundary 
values: A'{0) — and A'{oo) = 2. Using these values in Eq. (13), we obtain 
the gauge transformed flux $ = Anm/g. This does not imply, however, that 
the flux of the original vortex configuration $ is quantized in units of 'i'rr/g. As 
pointed out in Ref. [7], the singular gauge transformation is not really a gauge 
transformation at the origin. The term —d^ip on the R. H. S. of Eq. (14) injects 
a singular flux tube at r = 0. This is seen when one calculates the curl of this 
term as done in Eq. (12). The transformation of the gauge field given in Eq. 
(14) is accompanied by the transformation of the order parameter 

(j) ^ (/)' ^ cxp(—imip)(j) (15) 

Thus, whereas the vortex has a winding number to, the transformed vortex 
(j)' has a winding number 2to.. In both cases the flux remains to be quantized in 
units of 2Tr/g. 

The flux quantization for the vortex Ansatz (4) can be also deduced by 
relating $ to the topological (Pontryagin) index [6] 

g = -!- / e'^'n.idnri A d,.n)d^x (16) 

Stt J 

where n — ^^a$_ is the unit vector field representing the local orientation in 
the charge space. Applying Stoke's theorem to the R. H. S. of Eq. (8), we have 



$ - - / e^'^d^.A^d^x (17) 

With use of Eq. (4), we obtain n = (sin / cos mip, sin / sin rmp, cos /) yielding 

n.{df,n A dt,fi) = msin/(a^/9^^ - di,fd^,Lp) (18) 

From the definition of the gauge field A^ = —-^^9^^ we obtain using Eqs. 
(4) and (18) 

e^"9^i. = ^e^'^n.(a^nAa,n) - ^l + co&f)5\T^ (19) 

where the second term arises from the singularity of dip at r = (see Eq. 
(12)). Eq. (19) represents a generalization of the identity previously established 
for the case where zeros of the superfluid density (vortices) are absent [1],[8]. 
Using Eq. (19) in (17), wc have 

4> = -^ f e''''n.{d^n/\d^n)d^x = -—Q (20) 

Note that the second term of Eq. (19) docs not contribute to Eq. (20) since 
the boundary condition /(O) — n implies 1 + cos/(0) = 0. 

The quantity Q in Eq. (20) is the winding number which is the number of 
times the sphere S%^ij^ is traversed as we span the R^ space compactified into 



the sphere 5^,. It should be noted that this compactification is enabled by the 
boundary condition /(oo) = 0, implying 

lim n= (0,0,1) (21) 

r — >oo 

In this way, one is led to an i?^ space with infinity identified which is homo- 
topically equivalent to 5^ [6] . This equivalence is established by stereographic 
projection from the north pole of 5^^ to R^. The mappings Sgp -^ S'j^^u^ are rep- 
resentations of the second homotopy group 772(6*^) = Z where Z is the group of 
integers under addition [6]. Hence the degree of these mappings is Q = 0, ±1, ±2 
... The specific value of Q, consistent with the vortex Ansatz (4), is obtained by 
evaluating the local topological density in Eq. (16) from Eq. (18) in cylindrical 
coordinates 



1 Tn 

—e^^^n.id^nAd.n) = — sin fe^'" {d^fd,^ - d,.fd^^) 
OTT an 

m df m dcosf 

= - — sm/-- = -- ; (22) 

Anr dr Anr dr 

Using this result in Eq. (16), we obtain 

— dr — [cos/(oo) — cos/(0)] — —m (23) 

where we used the boundary conditions /(O) = tt and /(oo) — 0. Eqs. (20) 
and (23) imply $ = ^IL™ in agreement with Eq. (10). 

In conclusion, the present calculations show that the magnetic flux of the 
vortex in the SU{2) symmetric model of two-gap superconductor is quantized in 
the unit of — . For condensates formed by Cooper pairs, the coupling constant 
g can be related to the single electron charge e. We have g = || so that the 
unit — = -If = $0 where $0 is flux quantum carried by Abrikosov flux tubes 

[5]. 

The SU{2) rotational symmetry of the Lagrangian may not be easily achieved 
in real samples since it requires that the coupling constants of the potential term 
satisfy rather restrictive conditions [2, 3]. A more realistic potential energy has 
been considered in Ref. [1]. It has the generic form 

V = A + Bn3 + Cnl (24) 

where A, B and C are functions of masses and Copper pair (24) densities. 
The minimum of V takes place for n^ — n^ — —■^. Hence, the ground state 
value of the vector n is a circle on the 3%^^^ sphere. Now the SU{2) symmetry 
of the Lagrangian is broken to U{l)x U{1) symmetry. The two U{1) symmetries 
correspond to the conservation of integrals of the two quantities, A^i -I- N2 and 
Ni — N2, where Ni — |^ip is the density of the i—th condensate. 

If the potential in Eq. (1) is replaced by the U{1) symmetric form (24), the 
boundary condition (21) is replaced by 



lim n = I W 1 — rig cos (/3, w 1 — n^ sin ip, n^ ] (25) 



This implies that the infinity of the R^ space cannot be identified so that this 
space cannot be compactified into a sphere S'^p. Instead of the mapping S^p — > 
^"jieid^ we have a map of R^ to a fraction of the sphere 5*?^^;^ corresponding to 
the interval {n^ — —1,71,3). Consequently, the topological index Q becomes a 
fraction of to. Using cos/(oo) — n^ and cos/(0) = —1, Eq. (23) yields 

Q = -y(«3 + l) (26) 

so that the flux (20) becomes 

27r / 1 + n3 \ , , 

$.-to(^) (27) 

Consequently, the vortex in the U{1) model of the two-gap superconductor 
is quantized in units —a, where a = "^^^ is a continuous parameter satisfying 
the inequality < a < 1. 

According to Ref. [1] , the parameter B vanishes when the condensates satisfy 
the conditions toi = to,2 and A'^i = 7V2- In this case, the equilibrium value of n^ 
is fia = and Eq. (26) leads to a topological index Q — —m/2. Hence, we have a 
soliton which has half the winding number of the SU{2) vortex (mcron). At the 
infinity, the vector n in the mcron configuration points out radially in the xy- 
plane like the magnetization of a magnetic vortex in an easy plane ferromagnet 
[9]. We note that a soliton with Q — m/2 is also possible. Eq. (23) shows 
that this takes place when /(oo) = tt/2 and /(O) = 0. However, the boundary 
condition /(O) = requires a vortex Ansatz that is reasonable at the north pole 
of the S'?ieid sphere [10]. It is obtained by multiplying Eq. (4) by exp{im(p) 
thus reversing the winding number. Soliton with topological index Q = ±to./2 
is reminiscent of the meron excitations in the U{1) symmetric case in double- 
layer systems studied previously by Moon et al. [11]. In these systems the z- 
component of the pseudospin density at asymptotically large distances from the 
vortex center vanishes in order to minimize the charging energy of the capacitor 
formed by the two layers. This is formally analogous to our case where h = 
follows by minimizing the potential energy (24) where the parameter B is set 
equal to zero. However, the condition of symmetric bands (to,i — 7712, A^i ~ N2) 
does not seem satisfied in real materials. For instance, Ref. [12] confirms the 
existence of two different pair potentials and spatial scales for the two bands in 
MgB2 superconductor. 
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